Symmetry, Integrability and Geometry: Methods and Applications 



SIGIVIA 7 (2011), 111, 19 pages 



O 

(N 

o 

CD 

Q 
in 






in 

o 



Resolutions of Identity for Some Non- Hermit ian 
Hamiltonians. I. Exceptional Point 
in Continuous Spectrumtj 

Alexander A. ANDRIANOV^^ and Audrey V. SOKOLOV"^ 

' V.A. Fock Department of Theoretical Physics, Sankt- Petersburg State University, 
198504 St. Petersburg, Russia 
E-mail: andrianov@bo. infn. it , I avs^avs @rambler7rv\ 

^ ICCUB, Universitat de Barcelona, 08028 Barcelona, Spain 
E-mail: \andrianov@icc.ub. edu\ 

Received August 06, 2011, in final form November 25, 2011; Publislrcd online December 05, 2011 



littp://dx.doi.org/10.3842/SIGMA.2011.111 



Abstract. Resolutions of identity for certain non-Hermitian Hamiltonians constructed from 
biorthogonal sets of their eigen- and associated functions are given for the spectral problem 
defined on entire axis. Non-Hermitian Hamiltonians under consideration possess the contin- 
uous spectrum and the following peculiarities are investigated: (1) the case when there is an 
exceptional point of arbitrary multiplicity situated on a boundary of continuous spectrum; 
(2) the case when there is an exceptional point situated inside of continuous spectrum. The 
reductions of the derived resolutions of identity under narrowing of the classes of employed 
test functions are revealed. It is shown that in the case (1) some of associated functions 
included into the resolution of identity are normalizable and some of them may be not and 
in the case (2) the bounded associated function corresponding to the exceptional point does 
not belong to the physical state space. Spectral properties of a SUSY partner Hamiltonian 
for the Hamiltonian with an exceptional point are examined. 

Key words: non-Hermitian quantum mechanics; supcrsymmetry; exceptional points; resolu- 
tion of identity 
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1 Introduction 

The interest to exceptional points in non-Hermitian quantum dynamical systems has been re- 
voked recently [H [U [3l HI [5] although their very notion exists quite a time [6l O [HI [9l [10] . Their 
appearance can be associated to level coalescence at complex coupling constants for initially 
Hermitian Hamiltonians [11]. If existing they play an important role in definition of energy 
spectra and in construction of biorthogonal bases in Riesz spaces [12] . 

Whereas the appearance of exceptional points in discrete spectrum does not give rise to any 
principal obstacles for building a resolution of identity, the emergence of exceptional points inside 
or on the border of continuous spectrum makes the very construction of resolutions of identity 
rather sophisticated [2] . Their correct description in brief represents the main aim of our paper 
(entitled as Part I) whereas in the subsequent paper [13] (entitled as Part II) the detailed proofs 
of the results announced here are presented by one of us (A.V.S.). 

Let us start with the notion of exceptional point and further on outline the structure of the 
present work. The spectrum of a Hermitian Hamiltonian, in general, consists of continuous part 
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and discrete points. Meanwhile the spectrum of a non-Hermitian Hamiltonian may contain also 
a new type of spectral points embedded into a continuous spectrum and/or into a discrete one, 
namely, exceptional points. 

The exceptional point of the spectrum of one-dimensional Hamiltonian h defined on entire 
axis is an eigenvalue Aq of this Hamiltonian for which there is a normalizable eigenfunction 'ijjq{x) 
and also a number of associated functions i^jix), j = 1, . . . ,n — 1: 

ht/jo = XoTpo, {h - Ao)V'i = V'i-i, J = 1, • • • , ?^ - 1- 

For a discrete spectrum the latter ones are typically normalizable [21 [3] . On the other hand some 
non-normalizable associated functions bounded or even growing at infinity may be involved in 
building of resolution of identity as well. It will be proven in Part II. The number n is an 
algebraic multiplicity of Aq- Thus, in continuous spectrum one can deal with two types of 
algebraic multiplicities (which are not necessarily equal: their different types for continuous 
spectrum are discussed in conclusions). If an exceptional point Aq belongs to a discrete part of 
the spectrum then n simultaneously characterizes the order of a pole at E = Xq of the Green 
function. For an exceptional point Aq on the border of continuous spectrum the Green function 
reveals (see Section 12. 3p a branching point with the pole order 2n + 1 in the variable \/E — Aq, 
where n is a maximal number of linearly independent eigen- and formal associated functions 
of h for an eigenvalue Aq in the resolution of identity. When an exceptional point lies inside 
of continuous part of the spectrum the pole order may be larger than n (which has the same 
meaning as in the previous sentence) that is elucidated in details in conclusions. 

In this paper we build resolutions of identity for certain non-Hermitian Hamiltonians con- 
structed from biorthogonal sets of their eigen- and associated functions for the spectral problem 
defined on entire axis. Non-Hermitian Hamiltonians under consideration are taken with contin- 
uous spectrum and the following peculiarities are investigated: in Section [2] the case when there 
is an exceptional point of arbitrary multiplicity situated on a boundary of continuous spectrum; 
in Section [3] the case when there is an exceptional point inside of continuous spectrum. In Sec- 
tion [4] in conclusions the different ways to introduce algebraic multiplicities are discussed and 
the SUSY tools [Illll5l[l6l[I71lia[l9l[20l[2ll[22]in regulating them are inspected. 

More specifically in Sections [2] and [3] the reductions of the resolutions of identity under 
narrowing of the classes of employed test functions are elaborated. It is shown that the bounded 
associated function in an exceptional point inside of continuous spectrum does not belong to 
the physical state space (i.e. does not belong to the complete biorthogonal system built from 
eigenfunctions of the Hamiltonian and cannot be reproduced with the help of harmonic expansion 
generated by an appropriate resolution of identity). If an exceptional point lies on a boundary 
of continuous spectrum then some of associated functions included into the resolution of identity 
are normalizable and some of them may be not, still being elements of a rigged Hilbert space [23] 
and its dual one, the Gelfand triple generalized onto biorthogonal resolutions of identity. 

2 Resolutions of identity for model Hamiltonians 
with an exceptional point of arbitrary multiplicity 
at the bottom of continuous spectrum 

2.1 Basic constructions 

Let us consider the sequence of Hamiltonians, 

h^ = -d^ + '^^'^^ J xEM, d = —, Imz/0, n = 0,1,2,..., 

[x — zy ax 
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where /iq is the Hamiltonian of a free particle and all these Hamiltonians are 'PT-symmetric 
for the choice Rez = 0. One can easily check that for the Hamiltonian /i„ on the energy level 
E = there is an eigenfunction ?/^.„o(a^) and a chain of formal associated functions ^„;(rE): 

/in^nO = 0, hnl/jnl = 1pn,l^l, / = 1,2,3,..., 

{-^r(2n-l)n {-^r(2n-2l-m_ 

0!! = (-l)!! = l, {-2m-l)]\ = j^^^, m=l,2,3,.... (2.1) 

Moreover for odd n the functions il^niix); / = 0, . . . , [n/2] are normalizable (i.e. belong to L^(]R)) 
and when / > [n/2] they are non-normalizable and unboundedly growing for x — )• ±00. For even 
n the functions ipni{x); 1 = 0,..., [n/2] — 1 are normalizable, the function 

, , (-i)"(n-l)!! , , 

i;n,n/2ix) = ^ /^ „ (2-2) 

V^vrn!! 



is bounded but non-normalizable and the functions ipnlix) for / > n/2 are non-normalizable and 
unboundedly growing for x — )• ±00. 

The Hamiltonians hn, n = 0,1,2, .. . are intertwined by the operators 

qt=Td-Xn{x), xn(a;) = 7^ = ^, 9n V'nO = 0, n = 0, 1, 2, . . . , (2.3) 

with the help of the chain (ladder) construction [25l [261 123 [IZ| 

Kqt = Qthn-i, qnK = hn-iq~ 

and 

K = qnQni K-i = qnqn^ n = 0,1,2,..., h_i = -d'^ = Hq. 

One easily check that 

q~ipnl = -iipn-i,l-i, n,/ = 1,2,3, ... . (2.4) 

As well the eigenfunctions ^nix', k) for continuous spectrum of the Hamiltonian /i„ can be 
produced from the eigenfunctions 

'il;o(x;k) = ^=e''''', keW 
V2^ 



for the continuous spectrum of the Hamiltonian /iq of a free particle with the help of intertwining 
operators 



Mx;k) = -^l^l^%t-qte''', 

hn'ipn{x;k) = k^'ipn{x;k), /c G (-00, 0) H (0, +oo), n = 0,1,2, (2.5) 

For the function ^nix', k) one can derive the following explicit representation by induction, 

Mx; k) = ^ e^'^ ± ^V^,, /" ^ ^ (2-6) 

V2^ ^^ 2™m!(n-m !/c'^(x-z)'"' ^ ^ 
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where from it follows that the cofactor {i/k)"^ in (|2.5p provides the asymptotic form for 



'^n{x]k) 



ikx 



i + o( - 

X 



X — ;• ±00. 



(2.7) 



As well, using (|2.3|) and (|2.5|) . we can get the following representations for ijjnix; k): 



1 



Akz 



d n 



tpn{x;k) 



At last, it follows from (1231) and ^M that 



2tt \dt t 

ikz / d n 



d n — 1 
IE 



t 

d_ 
dk 



d 



1 



dt t^ 
n — 1\ / d 



k 



t=ik{x—z) 
1 



dk k 



ik{x—z) 



(2. 



'il^n{x;k) = -q^il)n-i{x]k) 
k 



Akz 



i{x — z) \dk k 



(^-l)['"^''^"-^(^=")]' ^ = 1'2>3,.... 



(2.9) 



Let us find now the connection between ij^nix; k) and ipniix), I = 0,1,2, ... . For this purpose 
we calculate the following derivative, using Leibnitz formula and the relation 4.2.7.14 from [28], 



am 

k^o dk''' 



[e-*^^A;>„(x;A:)] 



gn 



dk^ 



1 



M^-^) y" (n + g)! z^fc"- 



s=0 



2^s!(n-s)! (x-z)' 



fc=0 



27r2"(3;-z)"- 



inin{m,n} 



An+m 



j=0 



2tt {x - zY- 



IIP' 



m — 1) 



(-l)"m!^„_^/2(x), 

0, 

i™-"(m-l)!! 



, ^/2^(m - 2n - 1)!! (x - z)*^"" 
where C™ = n!/[r7i!(n — m)!] is a binomial coefficient. Thus 



77T, IS even, 

77T, is odd and ^ 2n — 1, 

?n is odd and > 2n — 1 , 



(2.10) 



V'ni(2;) 



1 ")" (9^' 

^ lim^[e-^'=^A:"'0n(a;;A:)]^ 



(2?)! k^odk'^^ 



1 = 0,1,2, 



(2.11) 



Since the continuous spectrum of the Hamiltonian /i„, n = 1, 2, 3, . . . coincide with [0, +00 
(see (|2.5p ). the eigenvalue £■ = of this Hamiltonian for the chain of functions ipni{x), I - 
0,1,2,... is situated at the bottom of continuous spectrum. It will be shown in Section 
which of these functions are included in the resolution of identity constructed from eigen- and 
associated functions of /i„. 

2.2 Biorthogonality relations 

The biorthogonality relations between functions ipniix), I = 0, . . . ,n — 1 follow from (|2.ip : 



tpniix)ipni'ix)dx = 0, l + l' f^n-1. 



(2.12) 



The biorthogonality relations between ipnoix) and ipn{x', k) can be derived with the help 
of (USD, (HID and ([2:9]) : 



'4'no{x) [k'^il)n{x; k)] dx = i i^noix) -d + 



+00 



n 



X — z 



[k^~'i;n-i{x;k)]dx 
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■- -i'tpno{x)[k'^ ^Tpn-i{x;k)]' 

n+OO 

+ i / ten V'n.o(a;)] [A;"" Vn-i(a;; A;)] dx = 0, n = 1, 2, 3, 



(2.13) 



The biorthogonality relations between normalizable associated functions ipni{x) and ipn{x:, k) 
can be derived in the same way with the help of (j2.4p by induction, 



ipni{x)[k""4'n{x; k)] dx = i il^niix) i -d -\ ) [P ^'ipn-iix;k)]dx 



oo 



X — z 



itpni{x)[k'' Vn-i(a;;A;)] +i [q^iljniix)][k'' Vn-i(a;; A:)] dx 

-oo J_^ 

X) 

ipn-i,i-i{x) [k"-~'^ipn-i{x; k)]dx = ■■■ 



-oo 

+ 00 



il^n-ifii^) [A:" ^il^n-iix; k)]dx = 0, 1 = 1,..., 



n — 1 



(2.14) 



The biorthogonality relations between non- normalizable formal associated functions ipni{x), 
I = [{n + l)/2], . . . ,n — 1 and ipn{x; k), 



ipniix) [k'^ipnix; k)] dx = 0, I 



n+l 



,...,n-l, 



(2.15) 



can be derived with the help of (j2.ip . (j2.6p . the Jordan lemma and the relation 4.2.7.17 from [28] 
as follows, 

X) 

ipniix)[k"-ipnix; k)] dx 



2l-n 



H) 



„(2n-2Z-l)!! t^ i"^(n + m)! 



E 



27r(2/)!! ^p2™m!(n-m)! 



+ 00 



(271-2Z-1)!! ,i^ {""{11 + my. 

^^~^' 27r(2n!! ^ 2™m!(n-m)!' 

. ..n (2n-2/-l)!! '^" J-(n + m)!_ 
^ ' (2l)\\ ^2^m\in-m)\ 

^ ' m=0 ^ ' 

n(2n-2/-l) 



{x 


- ^)2«-n-"^eifex ^^ 


oo 




71— m 


.+00 gifcx ^^ 


J 


_oo (a; - z)"-2'+m 


. d 
'dk 


\ 2l—n—m 


5{k) 
J 



+ sign (lmz)0(sign (Imz) k){—i) 
n 



Akz 



(2/)!! 



m=2l~n+l 



2™m!(n - m)!(n - 2? + m - 1)! 



2i-n 



k' 



{2n-2l-l)\\ ^^' J-(n + m)!_ 
^ ' (20!! ^2"^m!(n-m)! 

^ ^ m=0 ^ ^ 



"d^-^ 



(n— m) — 2(n— /) 



m 



+ sign (lmz)0(sign (Imz) A;)(— 1) 2 



„ (2 n-2/-l)!!(2/ + l)!! 2n-2^_. ,,, 
(2n)!! 



where 



^(t) 



m=2l—n+l 

1, O 0, 
0, t < 0. 
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The formal associated functions ipnlix), I = n, n + l,n + 2, . . . are not contained in the 
resolutions of identity (see Section [XBj) . but it is interesting that one can write the biorthogonality 
relations for these functions with ^„(x;A;) as well, 

) 

(-l)'i" -ikz\-^ i'"(n + m)! ''+°° 



27r(2/)!!(2/-2n-l)!! 



f-2-m!(n-m)! y„^ ^ 



I Z^ Om^\ n-mM \ dk ^ ^ 



(21)11(21 - 2n - ly.l ^^ 2"'m\(n - m)\ V dk 

\ / \ I m=0 ^ 

(20!!(2/-2n-l)!! ^p2-m!(n-m)! V^A: y* ^ ^'^ 

n 



(-1)" 



1)!! ^„ 2"*m!(n-m)! ^ '^ 



(2n!!(2/ - 2n - 1)!! '^ 2'"m!(n - m)! 

\ / \ I 771=0 



1 _^{2/-277)/^N y^ (^ + "i)K2^-"'-"i)' 



(20!!(2/-2n- l)!!(2/-2n)! ^ 2™m!(n-m)! 



777=0 



^-^,5(2^-2") (A;), / = n,n + l,n + 2,..., (2.16) 



where we have used (12.11). (12.61) and the relation 



E(n + m)!(s — m)! ^^ (n — TTi + 2m) {n — 1 + ■m)\{s — m)\ 
2'"m!(n-m)! ~ ^ 2'^m\{n - m)\ 

777=0 ^ ^ r77=0 ^ ' 

E{n — 1 + m)\{s — m)\ sr-^ {n — 1 + m)\{s — m)\ 
2™m!(n- 1 -m)! ^ 2'^-i(m- l)!(n-m)! 

m=0 m=l ^ ' ^ ' 

E(n — 1 + 7Ti)!(s — m)! ^^ (n + m)!(s — 1 — Tri)! 
2'"m!(n-l -m)! ^ 2'"m!(n - 1-m)! 

771=0 ^ '^ m=0 ^ '^ 

, v-^ (n — 1 + 77l)!(s — 1 — TTi)! , , , 

= (s + n > ^^ — ^-^ — ^ = --- = {s + n)---(s-n + 2){s-n)\ 

^ ^ Z^ 2™7n!(n-l-m)! v ; v )\ ) 

777=0 ^ ' 

= (s — n — l)!!(s + n)!!, s ^ n. 

At last, the biorthogonality relations between eigenfunctions for continuous spectrum of the 
Hamiltonian hn are proved in [T3] and take the following form: 

X) 

[k-^M^; k)] [(A:')"V'n(x; -A;')] dx = {k')'^''5{k - k'). (2.17) 

Let us notice that (|2J3]) - (fZT5|) contain ([232]) for / = 0, . . . , n-1, /' = and ^^TTJ) contains ([232]) 
for I = I' = due to (12. IIP . The relations ([2.16P can be derived with the help of differentiation 
from ([2T7]) also in view of ([ZTT]) . 

2.3 Resolutions of identity 

The initial resolution of identity constructed from ^„(x; A;) holds |13| . 

6{x — x') = / ipn{x;k)ipn{x';—k) dk, (2.18) 
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where £ is an integration path in complex k plane, obtained from the real axis by its deformation 
near the point A; = upwards or downwards (the direction of this deformation is of no difference 
since the residue of the integrand for the point A; = is equal zero in view of (j2.6p and (j2.10p ) 
and the direction of C is specified from — oo to +oo. This resolution of identity is valid for test 
functions belonging to CL^ = C|° n -L2(M; (1 + |a;|)'^), 7 > — 1 as well as for some bounded and 
even slowly increasing test functions (more details are presented in [13]) and, in particular, for 
eigenfunctions ipn{x',k) and for the associated function ()2.2p . 

One can rearrange the resolution of identity (j2.18|) for any e > to the forms 

5{x — x') 



') = 


( 1 — ^ /■+oo^ 

\J—oo Je / 

^X' - zV ^n-l^ 


j V„(x;fc)i/'„(x'; 
-.i{x;k)'il;n-l{x'; 


-k) 
-k) 


dk 

£ 


1=0 


yx-z) 


i{x — z) 




—e 



+ 



x' — 2;\"sinefx — x') 



X — z ) 7r(x — x') 



-e r+°°\ sme(x-x') 



+ / ) 'il^n{x;k)ipnix';-k)dk + 



7r(x — x') 



cose(x — x') Y^ (—1/4) 



n-l / . , .w mm{2l,n-l} 



'^ (-1/4)' """^^' '' {n+2l + l-m)\ ( x' - z \ ' 

i—n \ ' rr.—n ^ J \ / 



2'Ke{x — z)(x' — z) ^-^ e^Hx' — zY" 

^ '^ W=0 ^ ' m=0 

.ine{x-xr ^ {-l/4y -i^^-"-^> {n + 2l - m)\ ( x' - z ^ 



■kIx — z) ^-^ e^^ix' — z') 

^ > l=\ ^ > m=0 

^ m 

3=0 

n = 1, 2, 3, . . . , / = 1, . . . , 2n — 1, m = 0, . . . , min{/ — 1, n — 1} 

(cf. with (68) in [2]) and, consequently, to the form 



5{x - x') = lim'<. i + ] ipn{x;k)ipn{x';-k)dk + 



sin £(3; — x') 
7r(rE — x') 



_ cos£(x-xO "^ (-1/4)' "'°^^ '^ (n+2/ + l-m)! f x' - z ^ ™ 

27re{x-z){x' -z)f--e^^{x' -z^^ ^^ ^^+^''"'" (n - 1 - jn)! \x-z 

sinejx-x') ^ (-1/4)' --i^^'-^i^ (n + 2/-m)! f x^^z^) 

+ vr(x-z) A^e2/(3./_^)2/ Z^^ ^^^'-^'^^ri-l-m)l\x-zj J'^^'^"^ 

where the prime ' at the limit symbol emphasizes that this limit is regarded as a limit in the 
space of distributions (see details in [13] ) . 

From (j2.19p and (j2.20p one can try to derive various reduced resolutions of identity similar 
to the resolutions (69) and (70) of [2], which correspond to the partial case n = 1, or to the 
resolutions (15) and (16) of [29]. In particular, by virtue of Lemma 3.7 from [13], for test 
functions from CL-y, 7 > — 1 the resolution (|2.20p can be reduced to the form 

5{x-x') = lim'l( + \^n{x;k)il)n{x']-k)dk 

n-l , ^ , .w rmn{2l,n-l} 



cose{x-x') '^ (-1/4)' """^^' ^^^ {^n+2l + l-m)\ ( x' - z 

X — z 



27re(x — z){x' — z) ^ e^^{x' — z) 



^^ (-1/4)' """^^' ^V (^+2^ + 1-^)! A 

/=0 ^ ' m=0 ^ ' ^ 
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sins(x-xO :^^ (-1/4)' min{2Z_i.n-i}^ (^ ^ 2/ - m)! /x^\'- 1 ._^. 

+ 7r(x-z) A^e2/(2./_^)2/ ZL.^ ^2«,m,n^^_-^_^^, ^^ ^ _ ^ ^ J.^/./iJ 

More reduced resolutions of identity for the partial case n = 2 are presented in the forthcoming 
Section [231 

There is another way to transform the resolution (|2.18|) as well. This way was used for 
obtaining of the resolution of identity (26) from [29]. The integral from the right-hand part 
of ()2.18p is understood [13] as follows: 

iljn{x;k)iljn{x':,—k)dk= lim' / 'iljn{x;k)il)n{x';—k)dk, (2.22) 

where C{A) is a path in complex k plane, made of the segment [—A, A] by its deformation near 
the point A: = upwards or downwards (the direction of this deformation is of no difference as 
well as in the case with C) and the direction of C{A) is specified from —A to A. When using the 
facts that the integral in the right-hand part of (j2.22p is a standard integral (not a distribution) 
and that k'^""iljn{x;k)'ipri{x'',—k) is an entire function of k (see ()2.6p ) as well as employing the 
Leibniz formula and the formulae ()2.10p and ()2.1ip . we can transform ()2.22p as follows, 



il^nix; k)\ljn{x' ; —k) dk = lim' lim < I / + \ijjn{x]k)iljn{x;—k)dk 

A^+oo eiO y \J-A Je ) 



j=0 

' -e i-A 



lim' lim < / -|- / \i\)n(x]k)'\\)n{x\-k)dk 

A-!> + 00 £4-0 



n-l 



- 2(-l)" Y. (2n -21- lV2n.-2;-l Y. ^nm{x)^l^n,l-m{x) \ 

«=0 ^ ' m=0 ) 

= lim* < ( / +1 j 'il)n{x;k)'il^n{x';-k)dk 

n-l I "j 

- ^(-1)" E (2n -21- lV2n.-2^-l E i^nm{x)i^n,l-m{x') , (2.23) 

/=0 ^ ' ■m=0 ) 

where the latter equality is considered as a definition and the limit lim is regarded as point-wise 

£4,0 

one (not as a limit in a function space). 

Let us show that the terms outside the integral in the last line of (I2.23P can be rearranged 
as follows, 

n-l 



-2(-l)" Y ^2n -21- l)e^n^2l^l E i^n^{x)i^n,l-m{x') 
1=0 ^ ' m=0 

n-l 

= ^V'nKa;;e)V'n,n.-i-K3;';e); (2-24) 

1=0 

where ipni{x',s), I = 0, ... ,n — 1 is the chain of eigenfunction and associated functions (formal 
for I = [(n + l)/2], . . . , n — 1) of the Hamiltonian /i„ for the eigenvalue E = oi the form 

I 
ipniix;e) = 'Y'^ji^)'^n,i-jix), l = 0,...,n-l, (2.25) 

j=0 
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hnipno{x;e) =0, Kipniixie) = ipn,i--i{x;e), / = l,...,n-l 

with ttj (e), j = 0, . . . , n — 1 being unknown coefficients which will be found below. Using (|2.24|) 
and (|2.25p , it is easy to check that (|2.24|) is valid iff the coefficients aj (x) satisfy the following 
system, 

' 2f— !)"■ 

Y^a,{e)ai.,{e) = -j^^^j:^, / = 0,...,n-l. (2.26) 

After the redefinition 

«i(^) = ^2j fi ^J' J = 0, . . . , n - 1, 
where /3j, j = 0, . . . , n — 1 are new unknown coefficients, the system (|2.26|) takes the form 
^ 1 

^^^■^'-^■ = 2iTT' ^ = o,-..,n-i. 

The general solution of the latter system can be found in the recurrent form, 

1 1 / 1 '"^ \ 

The first terms of the sequence /3j, j = 0, . . . , n — 1 in the case /3o = 1 are the following ones, 

1 „ 31 ^ 863 „ 76813 



6 '^ 360' "^ 15120 '^ 1814400' 

Thus, we can choose the functions ^„;(2;;e) in the form, 

^nl{x\e) =i'^+^J-'^^Tpn,l-j{x), 1 = 0,. ..,71-1 

j=0 

and the resolution of identity holds. 



(2.27) 



5{x - x') = lim* < i + / ] 7pn{x;k)ipnix';-k)dk 

+ Y^ iini{x; e)V'n,„-i-z(x'; e) \ , (2.28) 

/=0 J 

(cf. with (69) from [2] for the case n = 1). Moreover, this resolution is equivalent to (j2.18p . i.e. it 
is valid for all test functions for which (12.18P is valid (cf. with the analogous results in Section [3] 
and in [29]). 

The resolution of identity (|2.28p contains all n functions from the chain ^„;(x;e), / = 
0, . . . , n — 1 and in this resolution for the eigenvalue E = there are no other eigen- or as- 
sociated functions of the Hamiltonian /i„. The order of the pole A: = for the Green function 



TTZ 

— V'n(a;>;A;)'0„(a;<;-A;) 
k 



'~ k=VE 

x> = max{x,x'}, x< = {x,x'}, Im\AE^O, {hn — E)Gn = S{x — x') (2.29) 
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considered as a function of fc = yE is equal to 2n + 1 in view of (j2.6p . and the exceptional 
point E = ^ coincides with the branch point of this Green function as a function of E. One can 
consider this pole of the order 2n + 1 as a result of confluence of the pole of the order n for the 
Green function as a function oi E = k'^ and of the factor k = \JE from the denominator of the 
Green function (see ()2.29p ). Thus, the number n of linearly independent eigen- and (formal) 
associated functions of the Hamiltonian /i„ incorporated in the resolution of identity (I2.28P for 
the eigenvalue E = ^ (exceptional point) is equal to the order of the "pole" i? = for the Green 
function Gn{x,x';E) in the sense elucidated above or, more rigorously, the order of the pole of 
Gn{x, x';E) as a function of fe = ^/E is 2n + 1 expressed in terms of the number n. 

Let us notice that in view of (j2.27p the functions iprii{x] e), / = 0, . . . , n — 1 satisfy the same 
biorthogonality relations from Section [2.21 as the functions ipni{x), / = 0, . . . ,n — 1. 

2.4 Example: case n = 2 

For the Hamiltonian 

6 

/i2 = -<9^ + 7 TT, X e M, Im z 7^ 0, 

{x — zY 

there are continuous spectrum eigenfunctions 

1 



ip2{x;k) 



1 



ik{x — z) k'^(x — z^ 



ikx 



h2ip2{x;k) = k'^4j2ix;k), A; G (-oo,0) U (0, +oo) (2.30) 



(see (j2.6p ) and also the normalizable eigenfunction ^20(2^) and the bounded associated func- 
tion ijj2iix) on the level E = 0, 



ll^20ix) = j=r- — , ^^21(2;) = -—7==, /l2'02O = 0, h2'4'21 = V'20, 



(see (|2.ip and (j2.2p ). In this case at the point k = there is a fifth order pole in the Green 
function G2{x,x';E) considered as a function of A: = yE (see (j2.29p and (j2.30p ) and thereby 
the exceptional point £■ = of the spectrum of the Hamiltonian /12 coincides with the branch 
point for the Green function as a function of E. 

In the case under consideration the biorthogonality relations (|2.12p - (|2.15p and ()2.17p take 
the form. 



-00 



00 

+00 



00 

+00 



ij^o{x)dx = 0, / ip2oix)[k^Mx;k)]dx = 0, (2.31) 

/+00 
V'2i(x) [k'^Mx] k)] dx = 0, (2.32) 

-00 

[k^Mx] k)] [{k'fi)2{x; -k')] dx = {k'f5{k - k'), (2.33) 

where (|2.3ip are included in (|2.33p due to the equality 

^2o(x) = lim rA:V2(2:;A:)l (2.34) 

fc— !>0 



(see ()2.1ip ) and (j2.32p can be derived from (j2.33p in view of (j2.34p and of the equality 

1 , 92 
2f-^o5F 



^2i(x) = -lim^^[e-^^^A;2v>2(x;fc)] 



(see (J21T]) as well). 
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It is straightforward to check that the resolution of identity ()2.19p (which is valid for test 
functions from CL^, 7 > — 1 as well as for some bounded and even slowly increasing test 
functions, see Section [2.3p . can be rewritten in the form 



5{x — x) 



— e /•+OON 

+ / "j ilj2{x;k)ip2{x';-k)dk 

00 J e 



r, / N, / / N , / ^ , / / NT sme(x-x') 6sin^ f (x - x') 
+ [i^2o{x;e)iij2i{x ■,e) + 'ip2iix;e)'ip2oix;e)\ H -, ^ + 



+ 



7r(x — x') ■K£{x — z){x' — z) 

12(x - x')sin^ ^{x- x')sin|(x- x') 3[e(x - x') - 2sin|(x- x')]^ 



7re^(x — z)2(x' — z)"^ 



+ 



2-K£'^{x — zY{x' — z)"^ 



(2.35) 



where the eigenfunction iIj2o{x]£) and the associated function ip2i{x]e) of the Hamiltonian /12 
read 

^20{X]£) = -i\ -i>20{x) = —=^ "2, 

V e V7re(x — zy 



V'2i(x;e) = -n/- 



^^21 (2;) + ^^'20(2;) 



2^^ 



1 + 



1 



e^{x — zy 



h2'il^2oix; e) = 0, /i2^2i(a;;e) =ip2oix;e). 



The eigenfunction "020(3;) s) and the associated function tlJ2i{x', e) obviously satisfy the biorthogo- 
nality relations similar to (j2.3ip and (|2.32p . 

It is shown in [13] that the resolution of identity (j2.35p can be reduced: a) for test functions 
from CL-y, 7 > — 1 to the form 



6{x — x') = lim'< 

e4,0 



00 J £ 



+00 



il)2{x;k)'4j2{x' ■,-k) dk 



6 sm \x — X 1 

+ [i)2o{x;£)i^2i{x';e) + '4)2i{x]e)'4}2o{x';£)\ ' ^ 



7re(x — z){x' — z) 



+ 



12(x - x')sin^ ^{x- x')sin|(x- x') 3[e(x - x') - 2sin|(x- x')]^ 



+ 



TTe^{x — zy{x' — zy 27re^(x — zy{x' — zy 

identical to (|2.2ip . b) for test functions from CL^, 7 > 1 to the form 



(2.36) 



6(x — x') = lim'< 



— e /■+OON 

+ / "j 'tp2{x;k)'ilj2{x';-k)dk 

00 J e 



+ [tp2o{x;£)tp2i{x';e) + ip2i{x;e)ip2oix';e)] 
12(x — x') sin^ ^{x — x') sin |(x — x') 



12(x — x') sin^ ^{x — x')sin|(x — x') 3[e(x — x') — 2sin|(x — x')]^ 



7re2(x — zy{x' — zy 27re^(x — zy{x' — zy 

and c) for test functions from CL^, 7 > 3 to the form 



(2.37) 



6(x — x') = lim'< 
£4-0 



+ / ^ ip2ix;k)ip2ix';-k)dk 

00 J e 



+ V'2o(a;; e)V'2i {x'; e) + ^21(2;; e)ip2o{x'; e) 



(2.38) 



12 



A. A. Andrianov and A.V. Sokolov 



The latter of these resolutions of identity seems to have a more natural form than the pre- 
vious ones, but the right-hand part of the latter resolution cannot reproduce the normalizable 
eigenfunction 



^P2o{x) CL^ = C^n L\R; (1 + |x|)^), 



7>3 



because of the biorthogonality relations. With the help of the Jordan lemma one can check that 



lim 



+°° 12(3; - x') sin^ 1(2; - x') sin |(a; - x') 



ip2o{x)dx 



lim 

£4-0 



7re2(x — z)^(x' — z)^ 



± le{z - xOe±^^(^-"') 



^^20 (2;') \ = -.i^-mix) 



(2.39) 



and 



lim 

£4-0 



+00 



3[e(x — x') — 2 sin |(x — x')] 



,/m2 



lim 



2'ke'^{x — z)2(x' — zY 
3 



--020(2;) dx 



±ie{z-x')l2 _^ ^_, _ ^r^±ie{z-x')/2 



^±ie{z-x') 



(2.40) 

'4'20{x') } = -'llj2oix'), 



where the upper (lower) signs correspond to the case Imz > (Imz < 0). Hence, just two 
last terms of the resolution of identity (j2.37|) and the corresponding terms in the resolutions of 
identity (j2.35p and (|2.36|) give a chance to reproduce ^^20 (2^) by these resolutions (cf. with the 
analogous results in Section 6.1 of [2], in Section 2 of [29] and in Section [3] of the present paper). 
It is interesting that contributions of these terms in the resolution of identity are (see Remark 3.4 
in |13j ) singular discontinuous functionals whose supports consist of the only element which is 
the infinity (cf . with the analogous comments in Section 2 of |29j and in Section [3] of the present 
paper). 

In the case under consideration the resolution of identity ()2.28p takes the following form, 



6(x — x') = lim* < 



+ / "j ip2{x;k)'ilj2{x';-k)dk 

00 J e 



+ ^2o(a:; e)V'2i {x; e) + ^21(2;; e)ip2o{x'; e) 



(cf. with (I2.38P ) and it is valid for test functions from CLj, 7 > — 1 as well as for some bounded 
and even slowly increasing test functions (see Section 12. 3p . 



3 Resolutions of identity for the model Hamiltonian 

with exceptional point inside of the continuous spectrum 



For the Hamiltonian 



2 nO^ix — z)sm2ax + 2cos'^ ax 

-o + 16a 



X G 



dx 



[sin 2ax + 2a(x — z)]'^ 
, a > 0, Iraz / 
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there are [2] continuous spectrum eigenfunctions 



ip{x;k) 



1 



1 + 



ik W'(x) 



W"{x) 



Akx 



A;2 - a2 W{x) 2{k^ - a"^) Wix) 
W{x) = sin2ax + 2a{x — z), 
hil){x; k) = k ip{x; k), A; G (— oo, —a) U (—a, a) U (a, +oo). 



(3.1) 



As well for the level E = a^ there is the normalizable eigenfunction ip^ix) and the bounded 
associated functiorlH V'i(2;)5 



V'o(x) 
-i/'i(x) 
such that 



(2a)3/2 



cos ax 



sin2ax + 2a{x — z) 
i 



^i{x) 



2a{x — z) sin ax + cos ax 



2x/2a 



g-*ax _ g*ax \^0\- 



2a [sin2ax + 2a{x — z)] 

X — )• iboo, 



(3.2) 



hipo = a'^tpo, (/i-a^)V'i = ^0- 

The exceptional point E^ = a^ is a pole for the Green function 



G{x,x';E) 



TTl 



— ijn{x>; k)ipn{x<; -k) 
k 



fc=v^' 



x> = max{x, x'}, x^ = {x, x'}, Im vE ^0, {h — E)G = 5{x — x'). 

This pole is a pole of second order, it is replicated on both sides of the cut -B > and there are 
no other poles for G{x,x';E). 

One can show [2] that the eigenfunctions and the associated function of h obey the biorthogo- 
nality relations, 



oo 

+ 00 



oo 

+ 00 



'il)^{x)dx = 0, / 'il)o{x)[{k'^ -a'^)tp{x]k)]dx = 0, 

/ + 00 
i^i{x) [(/c —a )ip{x; k)] dx = 0, 
-oo 

[{k^ - a^)tP{x; k)] [{{k'f - a^)i^{x; -k')] dx = [{k'f - a^f6{k - k') 
where (13. 3p are included in (j3.5p due to the equality 



'iI)q{x) = ^i\ — lim \{k — a )V'(2;; A:)] 
V a k^±a '-^ ' -' 

and (|3.4|) follow from (j3.5p in view of (j3.6p and of the equality 



(3.3) 
(3.4) 
(3.5) 

(3.6) 



^,(x) = ±H/-^hm|I-|-[(.^-a^)]^(.;.) 



1 =p 2iaz 
4a2 



V'o(x). 



The resolution of identity constructed from -(/'(x; k) holds [13], 
5{x — x') = I il){x;k)ip{x' ;—k)dk, 



(3.7) 



^There is a misprint in the normalization of ^o{x) and il)i{x) in [2]. 
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where L is an integration path in complex k plane, obtained from the real axis by its simultaneous 
deformation near the points k = —a and k = a upwards or downwards (the direction of this 
deformation is of no difference since for the points k = —a and k = a the sum of residues of the 
integrand is equal to zero). The direction of C is specified from — oo to +00. This resolution 
of identity is valid for test functions belonging to CL^ = C|° n L2(M;(1 + ja;!)'^), 7 > — 1 
as well as for some bounded and even slowly increasing test functions and, in particular, for 
eigenfunctions ^(x; k) and for the associated function ^1(3;). 

One can rearrange [13] the resolution of identity (|3.7p for any e £ (0, a) to the form 



5{x — x) 



—a—£ j-a—e r+oo 

00 J —a+£ J a+e 



tp{x; k)'ip{x'] —k) dk 



2 , sinefx — x' 

-\ — cosa(3; — X )- 



vr 



X — X' 



4a2 _ £2 

1 

vr 



cos 2a{x — X ) cos e{x — x) 



— '^q{x)iIjq{x') 
na 

2a 



l-2sin2-(x-x') 



4a2 _ £2 



sin 2a{x — x ) sin £(2; — x ) 



2"+e , dt 

['i/'o(x)'(/'i(2;') + ^i(x)^o(x')] / cost(x — x'^ 

2a-e 



t 



(3. 



and, consequently, to the form 



Six — x) = lim'< 



~a—£ j-a~£ j-+oo\ 

+ I + I ]ip{x]k)ip{x]-k)dk 



00 



a+£ J a+£ 

" 1 



2 , sme(x — x 

H — cos a[x — X ) 



vr 



X — X' 



e 



4a2 _ ^2 

1 

vr 



cos 2a(x — x') cos e{x — x') 



na 

2q 



1 — 2 sin — (x — x') 
2^ ' 



4a2 _ ^2 



sin 2a(x — x') sin e(x — x') 



[%l:o{x)i)i{x') + iPi{x)iIjq{x')] I cost(x-x') — 

2a-£ t 



(3.9) 



where the prime ' at the limit symbol emphasizes that this limit is regarded as a limit in the 
space of distributions. 

One can reduce [13] the resolutions of identity (|3.8p and (j3.9p for test functions from CL^, 
7 > — 1 to the form 



Six — x) = lim'< 
£4-0 



—a—e i-a—e r+oo 

00 J —a+e J a+e 



ipi^x; k)ipi^x'; —k) dk 



1 



vrea 



e 



1 — 2 sin -(x— x') iPq{x)iI)q{x' ) 



(3.10) 



and for test functions from CL^, 7 > 1 to a more simple form 



Six — x') = lim' < 
£4-0 



—a—e pa—e p+oa'\ -i | 

+ / +/ U{x-k)ij{x'--k)dk-—-^o{x)Mx')).{'i-ll) 

00 J — a+e J a+e / 



irea 



The latter of these resolutions seems to have a more natural form than the previous ones, 
but it cannot reproduce the normalizable eigenfunction 



Mx)^CL^ = C^nL\R;{l + \x\r), 7>1 
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because of the biorthogonality relations. With the help of (|3.2p . Lemma 4.8 from [13] and the 
Jordan lemma one can check that 



/ + 00 



■ sin^ -{x — x')iIjq{x)iI)o{x') 



irea 



Tpo{x) dx 



e4-0 



■ 2 ^ / /\ 
sm —{x — X 

■Kea 2 ^ ' 



2a 



cos^ ax sin2ax[sin2ax + 4a(a; — 2;)] 



(x — zy 



lim <^ ^0(3::') 



Aa'^{x — zf- 



■ 2 ^ / t\ 
sm —\x — X 

■Kea 2 ^ ' 



2a 



roi^) 



cos ax 

{x — z)"^ 



dx 



dx 



(3.12) 



lim 

£4-0 



e 



^ 



where the upper (lower) signs correspond to the case \va.z > (Imz < 0). Hence, just the term 

2 e 

svo? -{x - x')'iI)q{x)iPq{x') (3.13) 

Tiea I 



in the resolutions of identity (j3.8p - (|3.1Up gives an opportunity to reproduce 4'q{x) by these 
resolutions (cf. with the analogous results in Section 6.1 of [2], in Section 2 of [29] and in 
Section [2.41 of the present paper). It is interesting that the contribution of the term ()3.13p in the 
resolutions of identity ()3.8p - (|3.10p is a singular discontinuous functional (see Remark 4.2 in [13]) 
which support consists of the only element - the infinity (cf. with the analogous comments in 
Section 2 of [29] and in Section [2.41 of the present paper). 

There is another way to transform the resolution (|3.7p as well. This way was used for 
obtaining of the resolutions of identity (26) in |29] and (j2.28p in Section [2.31 The integral in the 
right-hand part of (j3.7p is understood as follows, 



ip{x;k)tlj{x';—k) dk = lim' / Tp{x]k)ip{x';—k)dk, 



(3.14) 



where C{A) is an integration path in complex k plane, obtained from the segment [—A, A] by 
its simultaneous deformation near the points k = —a and k = a upwards or downwards (the 
direction of this deformation is of no difference as well as in the case with C) and the direction 
of C{A) is specified from —A to A. Using 



(1) the fact that the integral in the right-hand part of ()3.14p is a standard integral (not 
a distribution); 

(2) the fact that (A;=p a)'^ip{x; k)ip{x'] —k) is a holomorphic function of fe in a neighborhood of 
k = ±a (see (IXTD l: 

(3) the Leibniz formula and the formulae (j3.2p : 

(4) the notation C{kQ\e) with fixed /cq G M and e > for the path in complex k plane 
defined by 

k = ko + e[cos(7r — {}) ±. isin(7r — '&)], ^ t9 ^ vr, 

where the upper (lower) sign corresponds to the case of upper (lower) deformations in C 
and the direction of C{ko;e) is specified from t? = to t? = vr; 
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we can transform (|3.14p as follows 
'ipi^x; k)tp{x'; —k) dk 



lim' lim 

A^+oo £4,0 







■ip(x; k)ip(x'; —k) dk 
dk 



j=0 

1 






i=o 



j\ dy 



a, 



^"4!{x;k)ij:n{x';-k)] 



r dk ) 



lim' lim 

A— !> + 00 £4-0 




—a—e j-a—e pA 

+ / +/ ]i'{x;k)ip{x;-k)dk 



+ 



1 



4:7ra 
1 



ipo{x)^Po{x') 



dk 



dk 



zvr 



lim' lim 

A-s>+oo e4,0 



C{-a;e) (^ + a)^ ^^(q;^) (^ " a)^ 

A; - a Jc{-a;e) ^ + « 



£(a;e) 




V'(x; k)'ip{x'; —k) dk 



1 



vrea 



^o(2;)'(/'o(a;') 




-a—e ra~s r+oo\ 

+ +/ ]7p{x;k)ip{x;-k)dk 

30 ^— a+e J a+e 



nea 



-M^)Mx')}, (3.15) 



where the latter equality is considered as a definition for lim* and the limit lim^io is regarded 
as a pointwise one (not as a limit in a function space). The resolution of identity (|3.15p is 
equivalent to (13. 7p . i.e. it is valid for all test functions for which ()3.7p is valid (cf. with (13. lip 
and with the similar results in Section 2 of [29] and in Section 12.31 of the present paper) . 

Let us note that the associated function tpi{x) does not appear in the derived resolutions of 
identity and is not expandable with the help of the resolution (j3.1ip . Thereby, this associated 
function does not belong to the physical state space (rigged Hilbert space). 

Let us notice also that the number (equal to 1) of linearly independent eigen- and (formal) 
associated functions of the Hamiltonian h for the eigenvalue E = a"^ included into the resolutions 
of identity p. lip and (j3.15p is less than the order (equal to 2) of the pole £■ = a^ of the Green 
function G{x, x'; E). 



4 Conclusions: indexes of exceptional points and SUSY 

We remark that, in general, one can introduce, at least, three different number indexes of 
exceptional point ii^ = Aq of a Hamiltonian h: 

(1) ni(Ao) to be a maximal number of linearly independent normalizable eigenfunctions and 
associated functions of h for the eigenvalue E = Xq; 

(2) 712 (Ao) to be a maximal number of linearly independent eigenfunctions and formal asso- 
ciated functions of h for the eigenvalue E = Xq appeared in the resolution of identity 
constructed from biorthogonal set of eigenfunctions and associated functions of h; 

(3) n3(Ao) to be an order of the pole in the point E = Xq for the Green function for h as 
a function of E (in the case with the Hamiltonian h = hn (see Section [2|), where the 
exceptional point E = Xq = coincides with the branch point of the Green function; it is 
natural to assume (see Section [2^3}) that ^3(0) = n). 
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It can be proven (by methods of [7j, see as well the example in Section 5.1 of [2J) that for an 
exceptional point outside of continuous spectrum all these indexes are identical, 

ni(Ao) = n2(Ao) = n3(Ao) 

and represent the algebraic multiplicity of the eigenvalue Aq. In the cases where an exceptional 
point is situated on the border of continuous spectrum or inside of it these indexes may be 
different. For example, in the case of an exceptional point E = at the bottom of continuous 
spectrum of the Hamiltonian /i„ in Section [21 



ni(0) 



n+1 



n2(0)=n, n3(0) = n => ni(0) ^ n2(0) = 713(0) 



and in the case of an exceptional point E = a^ inside of continuous spectrum of the Hamiltonian 
h in Section [3l 

ni(a ) = 1, n2(a ) = 1, 713(0 )= 2 =^ 7ii(a ) = 7i2(a ) < 713(0 ). 

Thus, one can consider these indexes as different generalizations of the notion of algebraic 
multiplicity. 

One can use SUSY technique [HI [151 US CZl 13 [13 [20] in order to regulate the algebraic 
multiplicity (in any sense mentioned above) of an exceptional point in the spectrum of a SUSY 
partner Hamiltonian with respect to the order of this exceptional point in the spectrum of a given 
Hamiltonian (originally proposed in [2T] and elaborated in details in [31[1]): 

(1) in order to increase the multiplicity of an exceptional point Aq in the spectrum of a SUSY 
partner Hamiltonian with respect to its order in the spectrum of a given Hamiltonian, 
one must take a formal eigenfunction (and a chain of formal associated functions) for 
the spectral value Aq of the latter Hamiltonian as transformation function(s) which tends 
(tend) to infinity for x — )• ±00; 

(2) in order to decrease the multiplicity of an exceptional point Aq in the spectrum of a SUSY 
partner Hamiltonian with respect to its multiplicity in the spectrum of a given Hamilto- 
nian, one must take a normalizable eigenfunction (and a chain of normalizable associated 
functions) of the latter Hamiltonian for eigenvalue Aq as transformation function(s). 

These statements can be clarified by the following simple example. For the eigenvalue 
i? = Aq = of the Hamiltonian /i„, n = 1,2,3,... in Section [21 there is a chain (|2.ip of 
the eigenfunction il)no{x) and associated functions ipniix), I = 1, . . . ,[{n— l)/2]. As well one can 
check that for the spectral value E = oi the Hamiltonian /i„, 7i = 0, 1, 2, . . . there is a chain of 
the formal eigenfunction (pno{x) and formal associated functions ipni{x), 

hn^nO = 0, hnfnl = 'Pn,l-l, / = 1,2,3,..., 

which tend to infinity for x — t- ±00. Thus, it can be easily found that a) if to use ipni{x), 
I = 0, ...,771 as transformation functions for /i„, then the resulting Hamiltonian is hn+m+i 
with the exceptional point E' = of larger algebraic multiplicity (except for the case with the 
indexes ni(0), tti = and odd 71, where the indexes 7ii(0) for /i„ and hn+m+i = hn+i are 
equal), and b) if to use ipni{x), I = 0, . . . ,m ^ [{^ ~ l)/2] as transformation functions for /i„, 
then the resulting Hamiltonian is hn-m-i with the exceptional point ii^ = of smaller algebraic 
multiplicity (except for the case with indexes 7ii(0), 771 = and even n, where the indexes 7ii(0) 
for hn and hn-m-i = ^n-i are equal). 
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